Bound states of a localized magnetic impurity in a superfluid of paired ultracold fermions 
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We consider a localized impurity atom that interacts with a cloud of fermions in the paired state. We develop 
an effective scattering length description of the interaction between an impurity and a fermionic atom using their 
vacuum scattering length. Treating the pairing of fermions at the mean-field level, we show that the impurity 
atom acts like a magnetic impurity in the condensed matter context, and leads to the formation of a pair of Shiba 
bound states inside the superconducting gap. In addition, the impurity atom can lead to the formation of deeply 
bound states below the Fermi sea. 



PACS numbers: 67.85.-d, 67.85.Lm, 37.10.Jk 

Magnetic impurities in superconductors are known not only 
to alter the BCS ground-state by introducing potential scat- 
tering, but also to be at the origin of the pair-breaking effect 
leading to elementary excitations fundamentally different than 
those present in pure superconductors. Their presence serves 
to attenuate superconductivity by formation of in-gap Shiba 
bound states 1 . In fact, at high concentrations magnetic im- 
purities induce gapless superconductivity 1 -. While magnetic 
impurities in a superconductor have often been discussed as 
one of the simplest models that exhibit interplay between su- 
perconductivity and magnetism, the consequences of this in- 
terplay are still not fully understood. 

The experimental realization of paired states in ultracold 
atom systems has shed new light on many problems in super- 
conductivity. In particular it has been instrumental in prompt- 
ing a better understanding of the BEC-BCS crossoveP^l in- 
clud ing t he case of pairing in systems with large spin imbal- 
anc d 15 l 16 l It has also prompted a new generation of research 
on the subject of the dynamics of theses system s 17 ! 18 !. 

Combine magnetic impurities and ultracold atom systems 
can have rich physical consequences, some of which we ex- 
plore in the present paper. In particular, introducing mag- 
netic impurities into fermionic superfluids would help in un- 
derstanding the interactions between magnetism and superflu- 
idity, and could help to resolve long standing problems such 
as how superconductivity becomes gapless. 

Alternatively, instead of studying how the system changes 
in response to magnetism, one can use localized magnetic 
impurities as a form of local probe. As an example, in the 
setting of high temperature superconductors, detection of the 
modulation of the local density of states by a magnetic impu- 
rity via a scanning tunneling microscopy (STM) was used to 
great advantage to probe the nature of quasi-particle states of 
these materials both in the superconducting and pseudo-gap 
phasePEl. 

Although STM spectroscopy is not currently possible 
in the ultracold atom setting Radio Frequency (RF) spec- 
troscopy^^ could be used to probe the nature of the su- 
perconducting state in the vicinity of the magnetic impu- 
rity. Combining a magnetic impurity with RF spectroscopy 
can be used to directly probe the size of the superconduct- 
ing gap eliminating the uncertainty due to effects like Hartree 



shift J 24 ! 2 ^ . Further, we envision that additional information 
from momentum resolved RF spectroscopy^ in the vicinity 
of a magnetic impurity could provide data on the symmetry of 
the gap and its nodal structure. 

In this paper we propose and investigate theoretically a 
scheme for introducing localized magnetic impurities into the 
ultracold atom fermionic superfluid. The impurity is formed 
by an atom of a different species (from the species making up 
the superfluid) that is localized by a deep optical lattice poten- 
tial. The laser frequency is chosen such that the optical lattice 
interacts only weakly with the two atomic species, |f) and \l), 
that make up the superfluid (i.e. |f) and \l) do not become lo- 
calized). The magnetic character of the impurity originates in 
the different interactions strengths between the impurity atom 
and |t) and \\.) atoms, which we describe by a pair of effective 
scattering lengths and a^. 

The main input into our theory of impurity localized states 
is the description of the atom scattering on a localized impu- 
rity. (1) We begin by showing that, under rather general con- 
ditions, we can describe the interaction between a localized 
impurity and the free atoms via an effective s-wave scatter- 
ing length. (2) As pointed out by Shiba, due to the sharpness 
of the BCS density of states, the magnetic impurity always re- 
sults in the formation of a pair of localized bound states, called 
the Shiba states. Indeed, we find that as long as ^ a^, the 
impurity atoms always induce a pair of bound states inside the 
superconducting gap. (3) Interestingly, we find that the Shiba 
state is not related to the under-sea bound state. By the under- 
sea bound state we mean the natural extension to the case of a 
filled Fermi sea of the Feshbach bound state formed between a 
fermion and a localized impurity in the absence of a Fermi sea 
when the effective scattering length is positive. Indeed, if the 
Feshbach bound state exists, it becomes the under-sea bound 
state when the Fermi sea is filled, remaining completely sep- 
arate of the Shiba state. (4) We show that both the under-sea 
bound states and the Shiba bound states can be resolved via 
RF spectroscopy. 

This paper is organized as follows: In section [I] we relate 
the bare scattering length between a pair of atoms in vacuum 
to the effective scattering length when one of the atoms is lo- 
calized by a parabolic confining potential. In section[ll]we use 
the effective scattering lengths to find the under-sea as well as 
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the Shiba bound states of the magnetic impurity. Next, we 
describe RF spectroscopy of the ultracold atom system with 



bound states in section III We discuss possible experimental 



realizations and atom species that could be used in section IV 



discuss the outlook in section[V] and draw conclusions in sec 
tionEl 



In order to relate the effective scattering length to the vac- 
uum scattering length, we must solve the scattering prob- 
lem. In general the scattering problem is complicated, and 
requires a numerical solution. In appendix [A] we state the 
scattering problem and derive an analytic solution for the spe- 
cial case of weak impurity-fermion interactions using a Born- 
Oppenheimer type approximation. 



I. EFECTIVE SCATTERING LENGTH 



The goal of this section is to show that the scattering of a 
fermionic atom off of a confined impurity can, under reason- 
able conditions, be described by a single quantity - the effec- 
tive scattering length. The problem of scattering on a confined 
impurity was previously studied in RefP^I here we review the 
basic arguments and summarize the results. 

We begin by assuming that the impurity-fermion scattering 
in vacuum can indeed be defined by a single scattering length 
for the s-wave scattering process. This condition means that 
the effective range ro of the impurity-fermion interaction po- 
tential is much smaller than the typical fermion wavelength 
1/kp, and thus we can treat r as being essentially zero. Since 
we want the fermion-impurity interaction to be tunable, we 
shall be primarily interested in operating in the vicinity of 
a wide Feshbach resonance (i.e. a resonance that meets the 
condition r <C 1/kp). If the effective range condition is not 
satisfied for the case of a free impurity (e.g. for the case of 
a narrow Feshbach resonance), it will not be satisfied for the 
case of a localized impurity, necessitating a more complicated 
description of the effective scattering process. Although, we 
do not treat the more complicated case in the present paper, we 
expect that the qualitative features, including the Shiba bound 
states, of a system with a narrow impurity resonance will be 
similar to those of a system with a wide resonance. 

In the problem with a confined impurity we have two im- 
portant energy scales: the typical kinetic energy of a scatter- 
ing fermion, which in our case is set by the Fermi energy scale 
ep, and the level spacing of the impurity atom which we label 
fajJi. We begin by pointing out that the scattering is elastic in 
the regime ep <C TtuJi. Further, in order for the scattering to 
be dominated by s-wave channel, we demand that the ground 
state wavefunction of the impurity must have a length scale 
yJKfm^uji that is much smaller than the wavelength of the 
scattering particle K] y/2m a ep (here to, stands for the mass of 
the impurity and m a for the mass of the scattering fermion). 
The two conditions are identical, up to a ratio of the masses. 
That is we demand that max(l, m a /mi)ep <C fiuii. 

Having derived the conditions for s-wave scattering we can 
write the resulting T-matrix for the scattering atom in the form 



T( W ) 



1 



(1) 



It is important to point out that since the impurity is local- 
ized the T-matrix features the fermion mass as opposed to 
the reduced mass /i = (m^ 1 + m^ 1 ) and an effective 
scattering length a a instead of the vacuum scattering ao, a . 



II. IMPURITY BOUND STATES 

In this section, we study the conditions for the existence of 
impurity bound states both in the normal (single component, 
non-interacting Fermi gas) and in the superconducting case. 
Our strategy is to obtain the T-matrix for scattering off of an 
isolated impurity in the presence of the Fermi-sea. Having the 
T-matrix, we can find the energies of the bound states from its 
poles. Further, we can also find the spectral function of the 
fermions, which we shall use in the next section to compute 
the RF spectra. 

In general, we can express the effect of the impurities on the 
Green function of the clean system G°(k, oS) via an expansion 
in the impurity density^ 

G(k,w) = G°(k,w) +n l G t) (k,w)T(u)G (k,uj) (2) 

+ 0(n?) 

where G(k, uj) is the Green function of the dirty system, T(u>) 
is the T-matrix, and rii is the impurity density. In this paper 
we shall always work in the dilute impurity limit, and thus 
drop terms of order 0(rij) and higher. The resulting equation 
is represented diagrammatically in Fig. [T^. T(uj) is obtained 
from the Lippmann-Schwinger equation 



T(u) = V + VQ Q {u)T(uS), 



(3) 



which relates the T-matrix to the impurity-fermion interaction 
potential V and the momentum integrated Green function of 
the clean system 



d 3 k 

(2^)3 



G°(k,w). 



(4) 



The Lippmann-Schwinger equation is illustrated diagrammat- 
ically in Fig.[lJ); it can be formally solved for the T-matrix by 
inversion 



T- 1 (u) = V- 1 -g°{cj). 



(5) 



Having specified the Lippmann-Schwinger equation, we 
first apply it to the case of an impurity in a one component 
non-interacting Fermi gas. This trivial case serves as an ex- 
ercise that demonstrates (1) regularization of point contact 
interactions, (2) properties of the T-matrix, and (3) relation 
between Feshbach molecules and under sea states. Having 
learned how to use the Lippmann-Schwinger equation in this 
context, we apply it to the T-matrix of the BCS state. 



3 



(a) 



X 



(b) 



FIG. 1. (a) diagrammatic representation of equation |2} and (b) of 
equation l[3j. Thin lines represent the clean (unperturbed) Green 
functions (G°(k,cj)), the thick lines the impurity-perturbed Green 
function (Gk), and the dashed line the interaction of these fermions 
with the impurity (V) 



The divergence in G°(u) is perfectly canceled by the renor- 
malized interaction to yield the T-matrix 



T( W ) 



m ( 1 
2tt \a 



1 

i^j2m{u 



(9) 



Unsurprisingly, the T-matrix has the same form as the vacuum 
T-matrix, but with frequency shifted by ep. This reflects the 
fact that energies must be measured with respect to the Fermi 
energy. The bound states of the system introduced by the pres- 
ence of the impurity are defined by the poles of the T-matrix. 
We thus find that a bound state exists only for positive values 
of the scattering length, with an energy 



uj b = -ep 



1 

2ma 



2 ' 



(10) 



A. Impurity in a one component non-interacting Fermi gas 



We first consider a fixed impurity interacting with a one 
component Fermi sea, the interaction being described by the 
scattering length a. Since the Fermi sea is non-interacting and 
the impurity is static, we can proceed simply by finding the 
one-particle eigenstates in the vicinity of the impurity poten- 
tial, and filling them up to the Fermi energy. For negative 
a, all eigenstates are part of the continuum, and there are no 
localized bound states on the impurity. As a becomes pos- 
itive, a single state, the Feshbach molecular state, with en- 
ergy — 1/ (2ma 2 ) peels off the continuum and becomes local- 
ized by the impurity (henceforth, the mass of the impurity no 
longer features and therefore we will use m for the mass of the 
fermion). When we fill the Fermi sea, the Feshbach molecular 
state appears as an under-sea bound state. 

In this subsection, we show how to recover this simple pic- 
ture in the T-matrix language. In the absence of impurity, the 
fermions are described by the following Green functional 



G(k,w) 



1 



u - 6c + iO+sgn(cj) ' 



(6) 



where £k = ek — £f = — £f and ep is the Fermi energy. 
In order to cancel the divergence of the integral of the Green 
function in the Lippmann-Schwinger equation we must use a 
renormalized interaction potential 2 -^! 



1 



2m 



V ATrh 2 a 



2m 



d 3 k 1 

(27r) 3 k2' 



(7) 



As described in Sec. [j] because the impurity atom is confined 
in the expression for the interaction potential we must use 
the fermion mass m and the effective scattering length o in- 
stead of the reduced mass and the vacuum scattering length. 
The momentum integrated Green function Q Q , that enters the 
Lippmann-Schwinger equation, can be obtained via contour 
integration 



0°M 



d 3 k 1 (2m) 3 / 2 , 

(2Try 6k 47T 



B. Impurity in BCS state 

In this subsection, we generalize the results of the previous 
subsection to the case of a localized impurity atom immersed 
in an ultracold BCS gas. We shall describe the BCS state 
at the mean-field level. Since BCS quasi-particles involve 
mixing particles and hole s, it is convenient to use Nambu's 
4-dimensional spinor basis^EIIIl 



*k = 



/ ' Ckt \ 

L -kt 



(id 



In this formalism, the BCS Hamiltonian becomes 



BCS 



6c -A 

o a A 

A 

-A -6c 



(12) 



where A is the BCS order parameter. The BCS Green function 
of the clean system is 



G°(k,w) = 



1 



W + 6sP3 + A(T 2 P2 



w - £,kP3 - Aer 2 P2 

uj + 6 
1 I W + & A 

k A n 



^ - il - a 2 

-A 




- & - A 2 I A u - £ k 











w - & 



(13) 



Here, {ai, a2, 03} and {pi, P3} are two sets of Pauli ma- 
trices, the first one operating on the spin space and the second 
on the particle-hole space. 

The interaction potentials between each of the two species 
that make up the BCS state and the impurity atom have the 
same form as the interaction potential in the single component 
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2m 



2m 



d 3 k 1 

(2tt) 3 k 2 



(14) 



Here corresponds to the effective scattering length between 
a |t) atom and the localized impurity, while between a |4_) 
atom and the impurity. In Nambu basis, the interaction poten- 
tial becomes 



V 



j_ 

V~2 



® P3 





1 

Vi 











\ 







i 
















1 

Vi 







I 











1 

v 2 


) 



(15) 



Substituting G°(k, lj) and V into the T-matrix equation |5]) 
we see that the four-dimensional Nambu space reduces into 
a pair of two-dimensional subspaces that can treated sepa- 
rately: the 'outer' (or 'first') subspace acts on the first and 
fourth Nambu components, whereas the 'inner' (or 'second') 
subspace acts on the second and third components. From this 
point, we will limit ourselves to one of them, say the first one. 

To use Lippmann-Schwinger equation in order to yield the 
T-matrix, our first step is the calculation of Q°, which can be 
written as the sum of a regular part Q® (lj) and a diverging part, 
as 



d 3 k 1 

From this definition, the regular part is 



(16) 



d 3 k 

(2tt) 3 



j 2 - - A 2 V ~ A w 



-A 



1 \ ( 1 

-1 



(17) 



Gj.(uj) can be expressed as a function of two integrals 



„ Lj 2 ~A 2 -( K 2 -e F ) 2 ' 
dn 



n 



LJ 



2 -A 2 -(K 2 -£ F ) 2 ' 



I 

where we have used the notation k = \z/y/2m. Both integrals 
can be evaluated using contour integration, to give 



(18) 
(19) 



h{v) = 
J 2 (w) = 



\J 1 e F + iVA 2 - lj 2 + yep - iy^A 2 



47rVA 2 - ( 

1 \Je F + i\/A 2 - lj 2 + y/e F - i\/A 2 



47rVA 2 - w 2 



+ A 2 -LJ 2 



(20) 
(21) 



Using the fact that that we take the correct branch of the square roots. Once 

more, the two diverging integrals in Q° and V^ 1 cancel, and 



^A^ 2 ~+Je F -i^A 2 ~^ 2 ~ fields 



1/2 



= 2 ^F + ve 2 F + A 2 -^ ^ (22) T _x M = ^ 2 ;: T j _ o (w) (24) 

we fj n( j Having solved the Lippmann-Schwinger equation, we can 

look at the properties of the resulting T-matrix. In particu- 

o /q / i b\ 1/2 

-n, , . w ' ( ep + c ) ' . . . . lar, we want to consider two regimes: bound states inside the 

y"lu}) = —t , -sgn(jR(a;)S(cLi)) , , , ., e 

rv ' 2-kJu 2 — A 2 gap and bound states outside the gap. 

c + ( £F -S) -A , (23) 



-A w-(eir-H) 
where S = \J e 2 F + A 2 — w 2 , and the sgn function ensures 
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Under-sea states 

Aiming to recover the Feshbach molecule-like bound state 
that we found to exist under the Fermi sea in the case of a one 
component gas, we make the approximation that A ~ 0. The 
bound state must correspond to a frequency uj = ujj, + i0 + , 
where ujb < — ep. Within this approximation, 



T-\u) 



o - 







2 iv a I 



2^ 



s/-e F - uj 








-iy/ep — uj 



(25) 



(26) 



We find that the T-matrix only has a pole (i.e. det T _1 (u>) = 
0) when the effective scattering length ay is positive. The fre- 
quency of the pole is 




FIG. 2. Representation of the sign of the solutions given by the first 
and second subspaces as a function of 1/kpa-f and l/hpa^. In the 
blue colored domain the solution given by the first subspace is nega- 
tive, while in the white colored domain the one given by the second 
subspace is negative. 



UJb 



-e F 



2ma 2 . 



(27) 



By looking at the complimentary 2x2 Nambu subspace, we 
find that another bound state exists for positive values of aj. 
with frequency uif, = — ep — l/2ma|. 

If we relax the approximation A ~ 0, we find that the 
under-sea state only becomes a sharp bound state in the limit 
ujb —> — oo. If the binding energy is not very large, then the 
under-sea bound state can serve as a Kondo impurity. How- 
ever, detailed analysis of this possibility is beyond the scope 
of the present article. 



Shiba states 




We now turn to the in-gap bound states predicted by Shiba, 
that is, | a; | < A. For weakly enough interacting BCS gases, 
we can make the approximation \u>\ < A <C e F , Within this 
approximation, 

/ m o \ 

\ U 27^ J 

+ WF=^ {-A uj J (28) 

The form of the T-matrix in the complimentary Nambu sub- 
space can be obtained from this one by making the substitu- 
tions A — » —A, ay <-> <Z|.. The poles of the T-matrix are 
defined by the equation 

U = ± 1 + kpaykp Oj, ^ 
VA 2 - uj 2 kpa^-kpa-f 

where the + sign corresponds to the first Nambu subspace 
and the — sign to the second Nambu subspace. From equa- 
tion ( |29l ), we see that as long as ay ^ there is exactly 
one pole of the T-matrix in each of the two subspaces. The 
two poles have opposite frequencies and correspond to the two 



FIG. 3. Energies of the two in-gap (Shiba) bound states as a function 
of l/kpa-f and l/kpa^ (here we use the approximation of \29) , we 
took A = 0.2ej?, and uj is measured in units of ep) 



Shiba states. We can interpret the negative frequency pole as 
a bound state for the quasi-particle of the gas, and the positive 
frequency solution as a bound quasi-hole. In Fig. [2] we split 
the {1/kp ay, 1/kpai} plane into two domains: the blue do- 
main corresponds to negative pole being in the first subspace, 
and the white domain to the negative pole in the second sub- 
space. The corresponding frequencies of the two Shiba states 
are plotted as a function of 1/kp ay and l/kpa± in Fig. [3] 

From ( p9| and Fig. [3] we see that the bound states are lo- 
cated inside the gap only for nonzero values of ay — 0|_. This 
fact can be quite straightforwardly interpreted: the interaction 
between the Cooper pairs and the impurity can be analyzed as 
the sum of a 'magnetic' term proportional to ay — a± and a 
non-magnetic term proportional to ay + aj,. The impurity can 
break Cooper pairs and give rise to in-gap states only when the 
magnetic term is finite. We note that when the non-magnetic 
term becomes zero, we recover the formula established by 
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l/kfa2=-0.5 



AA(fc,w) 



0.2 
U.I 
3 0.0 
-0.1 
-0.2 












— ^^fT" 








-4 -2 2 4 



l/kpa\ 

FIG. 4. Comparison between the approximate analytical solution 
(blue dashed curve) and the exact numerical solution (black curve) 
for the frequency (in units of £f) of the in-gap bound state (of the 
first Nambu subspace) as a function of 1/kpa-f, with A = 0.2eF 
and l/kpa_i = —0.5 



Shiba for a spin impurity in an electronic superconductor. 

Finally, we comment on the approximation that went into 
( |29| >. In figure[4]we compare the frequency of the Shiba state 
(of the first Nambu subspace) calculated using ([29]) and nu- 
merical solution of ( p4] i. We see excellent agreement between 
the approximate and exact answers, which persists to surpris- 
ingly large values of gap, A < 0.5eF- 



Discussion of bound states 

We underline that the Shiba and under-sea bound states 
are not related to each other. For example if both scattering 
lengths are negative, but unequal, then the two Shiba states are 
still present while the under-sea states are not. On the other 
hand for two positive and unequal scattering lengths there is 
a pair of under-sea bound states in addition to the two Shiba 
states. Finally if one scattering length is positive and the other 
is negative then there are again two Shiba states but only one 
under-sea state. 



III. RF SPECTROSCOPY 

We suggest that radio-frequency spectroscopy could be a 
good experimental probe for reading out properties of the 
Shiba as well as under-sea bound states. Basic tools for un- 
derstanding RF spectroscopy are given in 29 . RF spectroscopy 
works by converting |t) (or equivalently |t) ) atoms to a third 
hyperfine state labeled |3) by irradiating the system with pho- 
tons of frequency ui RF that bridges the energy difference be- 
tween |t) and 1 3) states. The bound states show up as edges in 
the spectra of transferred atoms when cjrf matches the bound 
state energy. 

In the following, we begin by reviewing the Fermi golden 
rule formula, in terms of |t) Green function, for the |t) — » 
1 3) transition rate as a function of wrf. Next, we apply the 




FIG. 5. Impurity induced correction to the spectral function 
A J 4(k, uS) of the one component Fermi gas for the case akp = 0.5 
(white - increase of spectral weight, blue - no change, red - de- 
crease), (a) Aj4(k, uu) as a function of momentum and frequency. 
The dashed white line indicates the position of the Fermi energy. 
Ayl(k, oS) shows a depletion of spectral weight along the clean dis- 
persion line k 2 /2m — £f (indicated by the red line), an under-sea 
bound state at ui — — 3ef, and excess spectral weight in the vicinity 
of the continuum band which corresponds to the impurity induced 
broadening, (b) Ayl(k, uj) as a function of frequency only with mo- 
mentum fixed at k — 0.5/cf [slice is indicated by the green line in 
(a)]. The spectral function can be decomposed into three (labeled) 
features: (la) a 5-function corresponding to the depletion of spec- 
tral weight along the clean dispersion line; (lb) part of the depleted 
weight is transferred into the vicinity of the clean dispersion line re- 
sulting in its broadening; (2) the remaining weight is transferred to 
a (^-function corresponding to the under sea bound state. We note 
that although the part labeled "Broadening" is divergent in the impu- 
rity density expansion, its frequency integral remains finite, and the 
spectral function fulfills the frequency sum rule. 



formula first to the case of one component gas and second to 
the BCS case. 



A. General formula for the RF transition rate 

We assume that the Hamiltonian of the system, subject to 
RF drive, may be written in the form 

H -^gas, impurity 

+ H 3 + Hkf, (30) 

where i? gas . impurity describes the fermion gas and the impurity, 
i?3 describes the Fermi ons in the |3) hyperfine state, and H R p 
describes the action of the RF radiation. In writing H in this 
form, we make the standard assumption that fermions in the 
|t) and |t) hyperfine states do not interact with fermions in the 
1 3) hyperfine state except through the action of H R p. Our goal 
is to calculate the RF current (i.e. the transfer rate of atoms 
from state |t) to state |3)) that is induced by H RF , which we 
do in second order perturbation theory (Fermi golden rule). 
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FIG. 6. Correction to the RF transition rate obtained for the one com- 
ponent gas due to the presence of impurities as a function of the drive 
frequency u, with k F a = —0.5 (top) and k F a = 0.5 (bottom). The 
RF spectrum for the clean case is sharply peaked at u — u-j, ~ with 
the width set by either trap properties and temperature. The impuri- 
ties have two main effects: (1) Since momentum is no longer a good 
quantum number, the impurities broaden the sharp absorbtion peak 
at uj — Co?, ~ 0. This broadening is composed of the depletion of 
the 5-function indicated by the blue arrow together with population 
of nearby-in-frequency states. (2) If there is a bound state, it induces 
an edge in the spectrum of transferred atoms followed by a broad 
feature indicated in pink. The broadening correction cannot be ac- 
curately captured in an expansion in impurity density. In fact at first 
order in impurity density we find that the correction is divergent but 
integrable. Therefore, in the figure we cut it off with a wavy line. 
While feature (1) is present independently of the sign of the scatter- 
ing length, feature (2) which corresponds to the coherent part of the 
transition rate correction (i.e. the bound state induced part) is present 
only for positive scattering length. 



The RF drive can be described by the Hamiltonian 



RF 



d 3 k 

(2tt) 3 



(^ RFt 4, kC T, k 



e^ 4 4 !kC3 , k ) 



(31) 

where il RF and w RF are the intensity and frequency of the RF 



drive; k (c3k) and cj k (c^k) are the creation (annihila- 
tion) operators for fermions in the |f ) and |3) hyperfine states. 
Since RF photons have a very small momentum (large wave- 
length) we neglect the momentum imparted on the atoms by 
the photons. Atoms in |3) hyperfine state are treated as free 
fermions and are described by the Hamiltonian 



H* = 



d 3 k 

(2tt)8 



(w 3 + e k )4 



3,k c 3,k, 



(32) 



where lu 3 is the splitting between the \f) and |3) states in vac- 
uum. The corresponding (Matsubara) Green function for |3) 



G 3 (k,iu n ) = - 



1 



iuj n - (e k + w 3 ) 



(33) 



The Golden rule formula states that current from |^) to |3) 



c 28 



/(cjrf) = 2fi 2 3 [D(iuj n -> w RF + i0+)] , 



(34) 



where 

V{iuj n ) 



/d 3 k 1 x 



(35) 

and uii and u> n are fermionic and bosonic Matsubara frequen- 
cies, respectively. Our Golden rule formula gives the tran- 
sition rate per unit volume. To obtain the transition rate per 
particle, we must divide /(wrf) by density [we shall use units 
where the density is set to k F /(6ir 2 ) = ^/2/(3n 2 )]. We re- 
state the golden rule formula in the more familiar real time 
version 

/<i 3 k de 
j^y^ A t (k, e)A a (k, e + w RF )n i ,(e), 

(36) 

where A a (k, cu) = — 23G (T (k, uj + i0 + ) are the spectral func- 
tions for a = {t, 3} fermions, (e) is the Fermi function for 
the f fermions, and we have assumed that the 3 band is empty. 
Using the fact that the |3) state is non-interacting, we can sim- 
plify this expression 



/(wrf) =n 2 



d 3 k 



3 A t [ k ' ^ + W3 - WRF 



x«f( n I-ws-wrf). (37) 

2 m 



Adding the assumptions that we are working at zero tempera- 
ture and the system has spherical symmetry, we can simplify 
the expression for the current even further 



y/2m(u^-^ 3 ) k 2 dk k 2 

_A t (k, — +u 3 -o JRF ). 

(38) 



j(wrf) = n 2 



To apply Eq. ( |36| l to the impurity problem, we separate the 
spectral function into that of the clean system A (k, to) and 
corrections that depend on the impurity density AA(k, uf) 



A Q (k,u) = A (k,u>) + 7i, [AA c (k,u) + AAi(k,u)} 



(39) 



Here, we have further separated the impurity contribution 
AA(k,w) = A A c (k,u)+AAi(k,Lj) into a coherent part that 
corresponds to the spectral weight of impurity bound states 
and incoherent part that corresponds to the broadening of the 
continuum states by impurity scattering. We apply the same 
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criteria to separate the RF transition rate 

I(u) = I (u) + rn(AI c (u) + AJ^cj)), (40) 

where I (ui) corresponds to the transition rate of a clean sys- 
tem, while AI c (uj) and AIi(u>) are the coherent and incoher- 
ent corrections due to the impurities. 

B. RF spectrum of a one-component gas with an impurity 

Suppose that the atom cloud is composed of a single, non- 
interacting, fermionic species in the hyperfine state |^) . To un- 
derstand the RF induced transition rate, and how it is affected 
by an impurity, it is useful to begin by describing the spec- 
tral function of the |t) fermions. The clean spectral function 
has the form A (k,uj) — 2ir5(ui — k 2 /2m + ep). The im- 
purity induced corrections to this spectral function AA(k, uj) 
are plotted in Fig. [5^. These corrections move spectral weight 
away from the clean dispersion and can be separated into an 
incoherent part that corresponds to the broadening of the con- 
tinuum band by impurity scattering and a coherent part that 
corresponds to the impurity bound states. 

In Fig. |5j3, we plot a slice through AA(k,uj) at fixed 
k = 0.5kp. In the slice we see three main features. First, 
we see a negative J-function feature, the location of which co- 
incides with the positive <5-function in Ao(k, 10) (feature la). 
This feature corresponds to the depletion of spectral weight 
from A (k, u). The spectral weight is transferred in to two 
regions: the under-sea bound state, which appears as a pos- 
itive 5-function in AA(k,u) (feature 2); the spectral weight 
is also transferred to the vicinity of the negative <5-function 
feature and corresponds to the broadening of the sharp disper- 
sion of the clean system (feature lb). Within our classification 
system, features la and lb correspond to incoherent spec- 
tral weight, while feature 2 corresponds to coherent spectral 
weight. Finally, we point out that although feature 2 is diver- 
gent, its frequency integral is finite. Indeed, the full spectral 
function satisfies the frequency sum rule, which means that 
the corrections satisfies 



= 



duj 
2tt 



(41) 



for all k. Having sorted out the spectral function we move on 
to the question of transition rate. 

Since the dispersions of the |f) hyperfine state and |3) state 
match, the clean part of the transition rate is sharply peaked at 
uj = 0J3 + ep and has the form 



on 



e F ] - cjrf). 



(42) 



At this point we pause to make several remarks. First, we re- 
mark that we have been following the notation in which the 
bottom of the |^) band is shifted to the frequency — ep. There- 
fore, the frequency difference between the bottom of the |f ) 
band and the bottom of the |3) band is lu 3 = lu 3 + ep. As a 
result, the frequency H13 and not uj 3 features in the transition 



rate formula Eq. ( |42| ). In the ultracold atom context, it is nat- 
ural to fix the "bare" splitting as Q3 instead of W3, since the 
bottom of the |f ) band does not move as the atom density is 
changed. Our second remark concerns the trapping potential. 
It is important to focus the RF radiation on the center of the 
trap in order to avoid the spatial smearing (due to shift of the 
Fermi energy), as discussed in RefP3. 

Next, we come back to the effects of the impurity. For pos- 
itive scattering length, there is an impurity bound state which 
results in a coherent correction to the transition rate 



A/ c (wrf) 



[uj 3 + e F }) - 1 



ma 2 (cjrf - [w 3 + e F )]) 2 



(43) 



In addition to this coherent correction there is also an incoher- 
ent correction, that occurs regardless of the sign of the scatter- 
ing length, and results in the broadening of the sharp transition 
rate of the clean state. We plot the impurity induced correc- 
tions to the transition rate in Fig.|6]for both negative and pos- 
itive scattering length. For the positive scattering length case, 
we highlight the coherent part of the transition rate, given by 
Eq. ( |4"3"j ), with pink shading. The incoherent part of the tran- 
sition rate correction is composed of a negative 5-function 
feature (indicated by an arrow in Fig. |6]l and a broad posi- 
tive feature. The (5-function feature corresponds to feature la 
discussed above: the depletion of spectral weight (and thus 
transition rate) from the clean spectral function. On the other 
hand the broad positive feature corresponds to feature lb: the 
broadening of the dispersion curve of the clean system. Since 
feature lb is divergent, we cut it off with a wavy line. As dis- 
cussed above, this divergence is a spurious consequence of the 
expansion in impurity density, and we do not expect to see it 
in experiment. 



C. RF spectrum of a BCS gas with an impurity 

In the clean BCS system, the fermion spectral function (for 
both species of fermions) has the form 



A),a(k,w) =— [(E k +i k )S(u-E k ) 



-(E k -Z k )S(uj + E k )}, (44) 



where E k = 



A 2 . The main feature of this spectral 
function is the superconducting gap in the density of states 
around the Fermi-surface. As before, the action of the impu- 
rity is to modify the clean Green functions and consequently 
the spectral functions. 

We begin by investigating how this spectral function is 
modified by the presence of the impurity atom, i.e. we com- 
pute -23G (k,w + iO+)T(uj + iO+)G {k,uj + i0+). We 
plot the change in the spectral function for both species of 
fermions induced by a magnetic impurity having kpa^ = 0.5 
and kpai = —0.5 in Fig. [7] Similar to the case of the sin- 
gle component gas, we see that the impurity has two effects. 
First, it induces a broadening of the continuum states. Second, 
it induces the formation of bound states. For the |f) fermions 



9 




FIG. 7. Impurity induced correction to the spectral function of the jf) atoms (left) and \\) atoms (right) as a function of momentum and 
frequency for the case a^kF = 0.5, a± = —0.5, and A/e = 0.4 (white - increase of spectral weight, blue - no change, red - decrease). The 
dashed white line indicates the position of the Fermi energy. Both A ni ^(k, uj) and J 4„ ii 4.(k, uj) show a depletion of spectral weight along the 
dispersion curve of the clean system indicated by the red line. ^4 nij -|-(k, uj) shows an under-sea bound state at uj ~ — 3e^ as well as a Shiba 
state at uj ~ — 0.13eF, while J 4„ ii 4.(k, uj) shows only a Shiba state at uj ~ 0.13eF- In addition, there is spectral weight in the vicinity of the 
dispersion curve of the clean system which corresponds to impurity induced broadening. 



it induces a Shiba state just under the Fermi energy, while for 
\\) fermions it induces a Shiba state just above the Fermi en- 
ergy. In addition, as is positive, the impurity induces an 
under-sea state for the |f) fermions that is analogous to the 
under-sea state of the one component gas. 

The RF spectrum for the clean BCS system is plotted in 
Fig. [8^, and the impurity induced corrections for the up and 
down atoms are plotted in Figs. [HJ? andJHJ;, respectively. The 
corrections to the RF spectrum due to the magnetic impurity 
are strongest for the |f) to |3) transition, depicted in Fig.|8]5. 
These consist of: (1) a dramatic filling of the gap, i.e. transi- 
tions to the left of the threshold frequency for the clean sys- 
tem, associated with the Shiba state below the Fermi energy; 
and (2) an edge in the spectrum that appears to the right of 
the main peak for the clean system associated with under-sea 
bound state. In the next three subsections we give analytical 
expressions for the RF spectrum of the clean system and the 
corrections due to under-sea and Shiba bound states. 



RF spectrum of the clean system 



associated with the band bottom. 



Under-sea states 



We follow the approximations of subsection II B oj < —ep, 
A«0, and use the approximate T-matrix of Eq. (|26]l. Around 
the pole uji, = —ep — 2j ^ a2 , the T-matrix takes the asymptotic 
form 



T(oj « uj d ) 



uj b 





o o 



(47) 



We recognize that in the vicinity of the bound state, the sin- 
gularity of the [1, 1] component of the T-matrix has the same 
form as the singularity of the T-matrix in the single component 
gas case, Eq. |9|. Thus, within our approximation AaO, the 
coherent part of the RF spectrum due to an under-sea bound 
state is identical to that of the single component gas, Eq. 43 
This contribution is indicated by the pink shaded region on the 
right of Fig. [HI). 



Using Eqs. ([36) and (j44j» we find that the transition rate for 
the clean system is 



/o(w) = : 



,m 3 / 2 AV(a)-w 3 ) 2 



2tt(w -ap- w 3 ) 5 / 2 



We note that by dividing our expression by the particle density 
we recover the transition rate per particle established by Ket- 
terle and ZwierleirP^. We plot this transition rate in Fig. [8^. 
The sharp onset at low frequencies corresponds to exceeding 
the threshold frequency 



OJth = w 3 



(46) 



Shiba states 



Following the assumption of subsection 



II B 



< A < 

(45) £p) and using the T-matrix of equation (|28|) we compute the 
coherent contribution to the RF transition rate from a Shiba 
bound state. For the coherent contribution we focus solely 
on the pole located at uj = — \bJb\, which exists in either the 
first or second subspace of the T-matrix depending on which 
domain of the {^-, ^} plane we are working, see Fig. |2j We 
assume that we are working at sufficiently low temperature 
so that only the negative frequency Shiba state is filled, and 
focus on the case of the negative frequency pole being in the 
first subspace. If it is in the second subspace, then the filled 
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Shiba state corresponds to a \\.) atom, and thus to detect it we Around the pole uib, the asymptotic form of the T-matrix is 
must use the RF transition \l) —> |3) instead of \f) — >• |3). found to be 



T(ui ~ uj b ) 



2?r 
mkp i i 



k F a-\ 



k F a 2 



1 

kpa2 



LJ — LOb 



-A 



-A 



LU — LUl, 



-R, 



(48) 



(49) 



where we define i? to be the regular part of the T-matrix in 
the vicinity of the pole. The coherent contribution to the spec- 
tral function must come from the above pole of the T-matrix. 
Combining the above form of the T-matrix with the clean BCS 
Green function Eq. 13 and the Golden Rule formula Eq. 43 we 
obtain 

AI c (u) = n 2mk ^ [G {k w u b ) -R-Goik^ut,)^ , (50) 

7T 



where k w — ^/2m(uj + ujb — ujs) and • indicates a matrix 
product and []n indicates the [1, 1] component of the matrix. 
From this expression, we see that for a Shiba state the thresh- 
old frequency for RF transition is cj t h = uj 3 — lj^- The coherent 
contribution of the Shiba state to the RF spectrum is indicated 
by the pink shaded region on the left of Fig. ^8}p. From the 
spectrum we see that most of the weight in the coherent part 
of the RF spectrum occurs at frequencies significantly higher 
than w t h- This is due to the fact that the Shiba state has most 
of its spectral weight concentrated at momenta ~ kp . 



IV. EXPERIMENTAL REALIZATION 

In this section, we turn to the experimental realization of 
such a system. In a typical dilute ultracold atomic gas, the 
Fermi wavevector will be on the order of kp ~ l/4000ao, 
where ao is the Bohr radius. The typical order of magnitude 
of the scattering length, in the absence of Feshbach resonance, 
is given by the Van der Waals interaction a ~ 50 — lOOao- In 
this regime, the kp(a^ — a^) and kp(a^ + a^) amplitudes 
always remain smaller than unity. Thus, in the absence of the 
resonance, the Fermion-Impurity (F-I) scattering lengths a+ 
and have roughly the same background values. As a result, 
the magnetic character of the interaction is vanishingly small 
and thus the Shiba states are too close to the gap edges to lead 
to observable results. 

The experimental conditions shall thus be chosen such as 
these two scattering lengths have widely different values, that 
is, close to an interspecies Feshbach resonance correspond- 
ing to one of the F-I interactions. Simultaneously, we wish 
to stay close above the F-F Feshbach resonance in order to 
maintain the large negative value of the associated scattering 
length. In conclusion, the impurity atom must be chosen to 



have a Feshbach resonance with one of the fermion hyperfine 
levels for a magnetic field slightly superior to the F-F reso- 
nant value. In addition to the requirement for Feshbach reso- 
nances, it is necessary to be able to confine the impurity very 
tightly in an optical lattice, while the fermions should still be 
relatively free. This would favor using a light fermion and a 
relatively heavy impurity atom, and employing a wavelength 
for the optical lattice that is near-detuned with respect to the 
optical transition of the impurity atom 

One possible choice of fermion atoms are the two lowest 
hyperfine states of 6 Li, which have a Feshbach resonance at 
Bq = 834 G. In order to achieve a BCS state, we want 
a "slightly superior" magnetic field, which means here that 
the difference between B and Bo shall be kept within the 
range of the Li-Li resonance width, which is approximately 
AB ~ 300 G. Amongst the few easily trappable bosons or 
fermions that could form a stable ultracold mixture with 6 Li, 
the boson 23 Na seems to rather well fit the above condition. 
Several Feshbach resonances have been observed between the 
23 Na hyperfine ground state and the ground state 1 1) of 6 Li, at 
magnetic fields close to the broad 6 Li- 6 Li Feshbach resonance 
[RefP^l. From these data, Gaesca, Pellegrini and Cote deduce 
in [Ref. 34 ] the existence of further resonances between Na and 
6 Li in states 1 1) and |2) between 834 and 1500 G. A co mplete 
list of predicted resonances is presented by Stan in Refs.EEH, 
along with a discussion of whether each corresponding hyper- 
fine mixture may or may not be stable towards losses due to 
spin-exchange collisions. For sodium-lithium mixtures, a pos- 
sibility is to use a green lattice laser at 532 nm. The effective 
mass for the sodium and lithium atoms as a function of lattice 
depth is plotted in Fig. [9^. For a lattice beam with ~ 120 /im 
waist, and a potential depth of about 4 lithium recoil energies, 
the sodium tunneling is essentially switched off (with an ef- 
fective mass of 77i* = 1000m), while lithium is still forming 
an itinerant Fermi sea (with an effective mass of m* ~ m). 

Another interesting combination are the lithium-rubidium 
interspecies resonances that have been found in RefpS. Here, 
there is a very interesting resonance at 882 G which is 1.3 G 
wide, not far from the 834 G resonance in lithium, and — 
as required for the assumptions in the paper — on the BCS 
side. The advantage of using the 882 G Rb-Li resonance over 
any of the Na-Li resonances is technical: the 882 G Rb-Li 
resonance has a width of 1.3 G while the Na-Li resonances 
have widths of ~ 300 mG. However, depending on the Li 
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FIG. 8. (a) RF transition rate for BCS state as a function of the 
drive frequency ui. Corrections to the transition rate for the \\) atoms 
(b) and |4-) atoms (c). (d) Total transition rate (clean+corrections) 
for 10% concentration of impurities, with divergences smoothed out. 
Throughout we have used a-fkF = 0.5, a± = —0.5, and A/e = 0.4. 
In (b) the coherent part of the transition rate correction, i.e. the part 
induced by the Shiba and the under-sea bound states is indicated by 
pink shading, with the peak on the left corresponding to the Shiba 
state and the peak on the right to the under-sea state. Similar to the 
case of the single component Fermi gas, the incoherent part of the 
transition rate correction is divergent at this order in impurity density 
(see Fig. [6J. Therefore, the total transition rate correction which is 
plotted in (b) and (c) is also divergent, and we cut it off with wavy 
lines, as before. 



density, the 882 G Rb-Li resonance may lie in the BEC-BCS 
crossover regime as opposed to the BCS regime. We suspect 
that the Shiba states will continue into the crossover regime, 
however determining their properties requires extending our 
theory. Alternatively, there is a Rb-Li resonance at 1067 G, 
which is wide (10.6 G), but lithium is then less strongly in- 
teracting, making it more difficult to attain a superfluid. For 
lithium-rubidium, one could use a laser tuned to about 820nm 
(see Fig. [9}}). As rubidium is so heavy compared to lithium, it 
makes for a very good localized impurity. 



V. OUTLOOK 

We suggest that the "implantation" of magnetic impurities 
into ultracold atom systems could lead to many exciting pos- 
sibilities. As already mentioned, one class of possibilities in- 
volves leveraging the interaction of magnetism and supercon- 
ductivity. This class includes the application of magnetic im- 
purities as local probes, which is the subject of the present 
paper. Another possibility is to study how the pair-breaking 
effect of the magnetic impurities leads to the destruction of 
superconductivity under various conditions. In 3D one would 
hope to realize the transition from gap-full to gapless super- 
conductivity. On the other hand, in 2D and ID the pair- 
breaking effect of the magnetic impurities is predicted to drive 
the superconductor-insulator transition. 

Another class of possibilities involves the Kondo effect. We 
already see a precursor to the Kondo level in the under-sea 
bound state. The nature of this under-sea bound state should 
undergo a dramatic transformation as we turn on the Kondo 
effect by changing the fermion-fermion interactions from at- 
tractive to repulsive. Significantly, using an optical lattice to 
localize the impurity atoms naturally invites the experimental 
realization of the Kondo-lattice model in the setting of ultra- 
cold atoms. The Kondo-lattice model is, in turn, a stepping 
stone on the path of studying itinerant magnetism. 

One significant difficulty in seeing magnetism in the setting 
of ultracold atoms has been the issue of achieving sufficiently 
low temperature. Perhaps magnetism without an underlying 
lattice could be technically advantageous. That is perhaps it 
will be easier to achieving the Kondo temperature by avoiding 
the lattice induced losses that feature prominently in the quest 
to achieve a magnetic transition (e.g. Neel temperature) in 
lattice systems. 



VI. CONCLUSIONS 

We have investigated the possibility of introducing a mag- 
netic impurity into a cloud of ultracold fermions. In particu- 
lar we have focused on the realization of a localized impurity 
atom that is immersed in a one or two component Fermi-gas. 
To understand the action of the impurity atom on the fermions, 
we have argued that it can be described by an effective scat- 
tering length, at least for the case of a broad resonance with a 
sufficiently tight impurity confining potential, which we relate 
to the impurity-fermion scattering length in vacuum. 

Using the effective scattering length description, we find 
the effects of the impurity on the free Fermi gas as well as a 
two component BCS condensate. In both cases we find that 
if there is a positive effective scattering length then the impu- 
rity forms an "under-sea" bound state. In addition impurity 
scattering breaks translational invariance and thereby broad- 
ens the spectral function of the clean system. Finally, for the 
BCS state if the impurity-fermion scattering lengths are dif- 
ferent then the impurity always induces a pair of Shiba bound 
states inside the gap of the superconductor. 

We demonstrate that the impurity bound states appear as ad- 
ditional features in Radio Frequency spectroscopy that should 
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FIG. 9. (a) Effective mass ra* of lithium and sodium atoms in a 532 nm lattice as a function of the lattice depth (measured in lithium recoil 
energies). Using a potential depth of ~ 4Er ] li it is possible to localize the sodium atoms (that serve as impurities) while lithium atoms 
remain itinerant, (b) Effective mass for lithium and rubidium atoms in a 820 nm lattice. Lattice depths between about 0.5 and 4 En^i can be 
used to localize the Rb atoms while Li remains itinerant. 



be detectable experimentally. Specifically, we suggest that 
6 Li BCS condensate with 23 Na impurities could be a poten- 
tially fruitful experimental system for study of magnetic im- 
purities. We speculate that beyond the study of bound states 
of dilute impurities, the same setup in combination with RF 
spectroscopy could be useful for study of gapless supercon- 
ductivity, Kondo effect, Kondo lattices, and other problems 
that combine localized moments and itinerant fermions. 
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Appendix A: Scattering Problem 

In this appendix, we state the scattering problem for the 
case of a confined impurity, and provide an analytic solution 
under special circumstances. We describe the scattering of a 
single fermion off of the trapped impurity by the Hamiltonian 



H 



Pi 
2m, 



Pn 



2m a 



(Al) 



where pi, and m, stand for the momentum, position, and 
mass of the impurity, p a , r a , and m a for the momentum, po- 
sition, and mass of the scattering fermion, Vi- a is a pseudo 
potential that describes scattering of the fermion off of the im- 
purity in vacuum. Eq. completely defines the scattering 
problem. However, in general, the equation must be solved 
numerically. 

For the sake of achieving a quantitate answer, we make ad- 
ditional assumptions. First, we assume that the effective range 
of the pseudo potential Vi- a is much narrower than the spa- 
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FIG. 10. Effective scattering length a that describes the scattering 
of a free fermion of mass m a on an impurity of mass m, localized 
in a harmonic potential of frequency Ui as a function of fermion - 
impurity atom interaction strength (scattering length in vacuum ao, a ) 
computed in the Frozen impurity approximation. For small ao,<*, a, 
depends linearly on ao, Q , Eq. l |A4fr , However, for large a we see a 
deviation from linear law. For large negative a it is possible to form 
bound states of the fermio n, which result in Feshbach resonances at 
(m a ao, a III) \J IrriiWi / it « { — 1.5, —9, ...}. 



tial extent of the harmonic oscillator ground state yjh/muii. 
Combining this assumption with the assumption that the typi- 
cal collision energy scales huji are much smaller than the char- 
acteristic resonance scale — % — , we can replace the interaction 



potential by a (5-function Vi~ a (ri — r a ) 



27rao ? , 



■S(ri - r a ). 



Finally, to obtain an analytical answer we make the frozen 
impurity orbital approximation. That is we first assume that 
we can write the wavefunction for the fermion and the impu- 
rity in a product form 



^(ri,r a ) = ipi(ri)ip a (r a ), 



(A2) 



and then we assume that i})i{ r i) is frozen to be the impurity 
ground state wavefunction. This approximation is similar in 
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spirit to the Born-Oppenheimer approximation that a heavy 
fermion is moving in the field of a fast impurity, with the 
additional assumption that the interaction between the two 
is sufficiently small that the impurity wavefunction is only 
weakly effected by the fermion. The approximation is valid 
for ao. a -C {n/mi)yJh/miUJi. The product wavefunction 
and the frozen impurity wavefunction approximations often 
appear in scattering theory. A particularly analogous problem 
where these approximations have been extensively used is the 
elastic scattering of a low energy electron from a hydrogen 
atorrP4 In this example, the role of the confinement potential 
is played by the electrostatic potential of the nucleus, which 
serves to localize the electron of the hydrogen atom. Within 
the frozen impurity wavefunction approximation, we define 
the effective potential that the fermion feels to be 



Vdr(r) 



(A3) 



where Tpo( r ) — (miUJi/Trh) 3 / 4: e~ muJir l 2Tl is the ground state 
wavefunction of the impurity. The effective scattering length, 
for small ao a> i s given by 



2.12 ao^ a m a 



1 + 



2 a , a m a y/miUi 



(A4) 



We note that even within our simple approximation, we find 
"geometric" resonances that are induced by bound states of 
V e g. These resonances can be clearly seen in the plot of the 
effective scattering length as a function of do a in Fig. 10 We 



expect that some of these resonance would survive in a more 
complete theory of the scattering process, and could be used 
in experiment to tune the "magnetism" of the impurity. 
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